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We investigate the effects of mesoscopic inhomogeneities on the metal-superconductor transition 
occurring in several two-dimensional electron systems. Specifically, as a model of systems with 
mesoscopic inhomogeneities, we consider a random-resistor network, which we solve both with an 
exact numerical approach and by the effective medium theory. We find that the width of the 
transition in these two-dimensional superconductors is mainly ruled by disorder rather than by 
fluctuations. We also find that "tail" features in resistivity curves of interfaces between LaAlOa or 
LaTiO:? and SrTi03 can arise from a bimodal distribution of mesoscopic local T c 's and/or substantial 
space correlations between the mesoscopic domains. 
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I. INTRODUCTION 

The recent discovery of superconductivity in transition 
metal oxide interfaces 1 - - — and of pseudogap effects and 
disordered-induced inhomogeneity in thin conventional 
superconducting films 4 ^— has triggered a renewed interest 
in transport and superconductivity in two-dimensional 
electronic systems. 8 Owing to their low dimensionality 
and to unavoidable defects in the growing processes, dis- 
order plays an important role in the physics of these 
systems. First of all, impurities lead to localization ef- 
fects, which can degrade the superconducting critical 
temperature^ and can give rise to localized preformed 
Cooper pairsJ^r— As it has become evident in a num- 
ber of recent numerical works ; 13 ' 14 the system manifests 
a spontaneous tendency towards inhomogeneities even in 
the presence of a homogenous disorder distribution. In 
particular disorder can also give rise to patchy electronic 
systems which are inhomogeneous on a mesoscopic scale, 
with superconducting regions embedded in, or connected 
by, normal/insulating regions ^ 15 ' 16 

Understanding all such issues is a difficult task, partic- 
ularly close to the superconductor-insulator transition, 
and it would obviously involve together concepts like 
percolation, establishment of coherence between super- 
conducting "islands" , role of disorder in the Berezinski- 
Kosterlitz-Thouless transition, and so on. In this con- 
text it would be useful to disentangle the various as- 
pects by understanding whether and which properties of 
these systems could be explained just in terms of meso- 
scopic inhomogeneities. Therefore we find it important 
and timely to investigate systematically the effects of 
large scale inhomogeneities on the superconducting tran- 
sition in two dimension, irrespectively of their micro- 
scopic origin. The analysis of these effects seems particu- 
larly compelling in metal-oxide interfaces, upon inspect- 
ing their sheet resistance Ra(T) data around the criti- 
cal temperature. While homogeneously disordered films 
of conventional superconductors present rather sharp 
transitions with high slopes of i?n(T) around T c ^~— 




FIG. 1. Sketch of the typical resistance curve for a supercon- 
ducting interface (circles). The width of the transition can be 
estimated as AT C /T C ~ (Th — Ti)/Ti. Notice the persistance 
of a pronounced tail also below T; . By closer inspection of the 
data in Ref. 0, one can see that for negative applied voltages 
< —80 V such a tail saturates to a low but finite resistance 
value. 



these superconducting heterostructures commonly dis- 
play broad transitions 2 ^ 3 - even for relatively small val- 
ues of the normal-state resistance. The typical shape 
of near the transition is sketched schematically 

in Fig. [1] In general, the downturn of R^(T) towards 
zero is characterized by a linear regime with a relatively 
small slope. By defining the two scales 1) and Th as in 
Fig. [TJ one can estimate approximately the width of the 
transition as AT/T C w (T h - T/)/Tj. The typical val- 
ues for AT/T C for superconducting interfaces are listed 
in Table I, along with the corresponding values in dis- 
ordered films of conventional superconductors, as TiN, 4 
and NbN. 7 Inspecting Table I one can easily see that the 
typical width of the superconducting transition in metal- 
oxide interfaces is substantially larger than the width of 
standard disordered films, that attain values compara- 
ble to those reported for the interfaces only at extremely 
large disorder concentration, where also the critical tern- 
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Ref. 


Sample 


T h (K) Ti (K) AT/T C 


[2] 


V=180 V 


0.347 


0.263 


0.32 




V=100 V 


0.338 


0.260 


0.3 




V=0 V 


0.306 


0.2 


0.53 




V=-40 V 


0.26 


0.138 


0.88 




V=-60 V 


0.235 


0.100 


1.35 


® 


TiNl 


1.71 


1.34 


0.28 




TiN2 


1.34 


1 


0.34 




TiN3 


0.83 


0.47 


0.76 


W 


k F l = 3.27 


8.27 


8.50 


0.03 




k F £ = 1.58 


2.36 


4 


0.69 



TABLE I. Estimate of the superconducting transition width 
in superconducting interfaces (Ref. 0) and disordered films 
of conventional superconductors (Refs. Q for TiN and for 
NbN). The labels used to identify the samples are the same 
used in the corresponding publications. In the case of Ref. 
V is the bias potential applied to the SrTiC>3 substrate to 
change the charge density at the interface. 



perature is almost driven to zero. The observed tempera- 
ture dependence of the sheet resistance in superconduct- 
ing interfaces shows instead two characteristic features 
not reported in disordered films: (i) the transition width 
becomes large for relatively small variations of the critical 
temperature (compare the case V = —40 V and V = —60 
V in Table I) ; (ii) an extended tail is observed below 7} in 
a sample with an applied bias V < — 80V, with the resis- 
tance saturating to a low but still finite value.— These two 
features call for a specific analysis to see whether in het- 
erostructures these large transition widths and persisting 
tails arise from mesoscopic inhomogeneities. 

It is worth noting that such broad transitions can 
hardly be ascribed to ordinary superconducting fluctu- 
ations. Indeed, considering the Aslamazov-Larkinii con- 
tribution of superconducting fluctuations to the paracon- 
ductivity [see Eq. ((5|) below], one obtains that AT/T C — 
Rn/Rc, where Rn is the normal-state sheet resistance 
(i.e., approximately the value right above the beginning 
of the downturn of i? n (T)) and R c = 16h/e 2 = 65.6 kfl. 
Since Rn values for the interfaces are of order of 1 kO 
or even smaller, ordinary superconducting fluctuations 
would give AT/T C ~ 0.02, i.e., at least one order of mag- 
nitude smaller than the one reported in Table I. This re- 
sult has to be contrasted, for example, to other systems 
like high-temperature superconducting oxydes^ where 
the conclusion was reached instead that fluctuation ef- 
fects are enough to account for the width and rounding 
effects around T c , while mesoscopic inhomogeneities play 
a minor role. 

The disordered mesoscopic model we have in mind con- 
sists of large metallic regions with randomly distributed 
critical temperatures that we map into a network of ran- 
dom resistors. Upon decreasing T more and more re- 
gions (links) become superconducting until a percolation 



threshold is reached and the superconducting transition 
occurs. The model can be solved by means of the so- 
called effective medium theory (EMT) , 19 ' 20 that is of- 
ten used as a valuable guide to investigate the effects of 
mesoscopic inhomogeneities. However, to establish the 
degree of reliability of EMT in the case of different disor- 
der realizations, we will compare the EMT results with 
exact numerical solutions on finite clusters. In particu- 
lar, owing to the "mean-field-like" character of EMT, it 
is important to establish to what extent this approach 
can reliably be used in cases where space correlations are 
sizable. 

We notice that in our model the superconducting re- 
gions are assumed to be large enough to have a fully 
established local coherence and to make negligible the 
charging effects. As a consequence, neighboring super- 
conducting resistances immediately establish a mutual 
phase coherence as soon as both have become supercon- 
ducting. This clearly distinguishes our framework from 
the case of granular superconductors, where the grains 
have usually nanoscopic sizes of the order of the coher- 
ence length of the pure system. Moreover, the large size 
of a mesoscopic domain also allows to consider paracon- 
ductivity effects (both d la Aslamazov-Larkin 17 and a la 
Halperin-Nelson^) to occur within each domain. 

Our paper is structured as follows. In Sec. II we sum- 
marize some results of the EMT in the case of mixed 
normal-superconducting states. In Sec. Ill we describe 
our numerical approach, where we find an exact solution 
to a network of random resistors with different local resis- 
tances and/or critical temperatures. Sec. IV is devoted 
to the analysis of different effects (T c distribution, inclu- 
sion of correlation between local T c and local resistance, 
gaussian or Berezinski-Kosterlitz-Thouless paraconduc- 
tivity fluctuations) in the absence of space correlations. 
This last limitation is overcome in Sect. V, where differ- 
ent realizations of space correlations are considered. Our 
concluding remarks are reported in Sec. VI. 



II. EFFECTIVE MEDIUM THEORY 

The effective medium resistivity p em of a random re- 
sistor network, in which the value pi of the resistivity of a 
resistor on a bond of the network occurs with a frequency 
Wi, obeys the equatio n 19 i 20 
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where the sum is carried over all the possible values of 
Pi, ^2iWi = 1, and the parameter a is related to the 
connectivity of the network. For a cubic network in D 
spatial dimensions, a = D — 1. Eq. (JXJ) can be recast in 
the self-consistent form 
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With p P < ${p em ) < PS, PP = [J2i W i/Pi\ ± i PS = 

^ i w i p i , and <&'(p em ) > 0. It is therefore evident that 
the function $(p e m) interpolates between the parallel 
pp and the series ps of the resistors. Since in Eq. (fT} 
-1/a < F{p em ) < 1, and T'{p em ) > 0, Eq. fll} has a 
unique solution, which can be efficiently found, e.g., by 
bisection, within the interval [p m in, Pmax], where p m in 
and pmax are, respectively, the smallest and largest pos- 
sible values of the resistivity. Indeed, from Eq. (J2j), it is 
evident that p mm < p P < p em < p s < p M AX- We spe- 
cialize the EMT to the case in which the random resistor 
network undergoes a metal-superconductor transition, by 
varying some control parameter, e.g., the temperature T, 
on which the resistivities pi depend. Lowering the tem- 
perature on the metallic side, an increasing number of re- 
sistors become superconducting and p em decreases. Near 
the transition p em — > + , and Eq. (fT]) can be linearized 
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whence it is evident that the transition occurs when 
the total weight of the superconducting resistors, w s = 
J2 Pi =o Wi > ecm als 1/(1 + a). On a square lattice a = 1 
and we find that the EMT correctly reproduces the per- 
colation threshold w s = 1/2 in two dimensions. 

To gain further insight into the physics of the metal- 
superconductor transition within the EMT, let us ini- 
tially assume that the resistivity on a bond of the ran- 
dom resistor network may take only two constant val- 
ues, and po > 0, and that the network is charac- 
terized by a distribution of T c 's W(T C ) of the individ- 
ual resistors, with /_ °° dT c W(T C ) = 1. Henceforth, 
to simplify the notation, we adopt a continuous de- 
scription of the statistical distributions. Each resistor 
has a resistivity po at high temperature and its resis- 
tivity vanishes as soon as the temperature T is low- 
ered below its critical temperature T c , i.e., the resistiv- 
ity of a bond is taken as p = poO(T — T c ), where po 
is the same for all bonds and T c is a random variable. 
The distribution of resistivity in the system is there- 
fore w(p) = w s {T)S{p) + [1 - w s {T)]S{p ~ po), where 
w s (T) = J^°° dT c W{T c ) is the statistical weight of the 
superconducting resistors at a temperature T, i.e., the 
frequency of occurrence of resistors with T c > T. In such 
a case, Eq. (fTJ) may be readily solved, to yield 20 



p em (T) = (1 + a)po / dT c W(T C ) (3) 

for T > T", and p em (T) = for T < T C Q , where T" is 
the critical temperature of the effective medium, which 
is defined by the equation 
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Thus, for T — > oo, p em — > Po, and p em (T c Q ) = 0. On pass- 
ing, we note the interesting inversion formula W(T C ) = 
p'em( T )\T=Tj[{l + &)po\, which holds for T > T? , and 
allows one to reconstruct the distribution of T c 's for 
T c > T", from the behavior of the effective medium re- 
sistivity as a function of T. This relation holds provided 
that the normal-state resistivity po does not depend on 
temperature. We apply the analytical solution ([3]) to a 
paradigmatic example which shall be used as a bench- 
mark in the forthcoming analysis. Let us assume a dis- 
tribution of TVs 
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with W\ + w 2 = 1. This is a bimodal distribution of T c 's, 
with two characteristic values, T\ and T 2 , and widths a\ 
and (72. As a limiting case, e.g., for w 2 — 0, a single 
Gaussian distribution is recovered. Specializing Eq. Q 
to the present case, we find 
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where erf(x) = f dze 
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As it will be clear below, such bimodal distribution turns 
out to reproduce quite well a large tail in the resistiv- 
ity as T" is approached, as in Fig. [TJ Indeed, while 
for a single gaussian distribution both the EMT and the 
random-resistor network solution give a resistivity van- 
ishing almost linearly at T", when part of the system 
is not superconducting (i.e., T 2 vanishes) the transition 
occurs smoothly with a resistance vanishing with an up- 
ward curvature, in close resemblance with the experimen- 
tal data in superconducting interfaces 



III. EXACT SOLUTION OF THE 
RANDOM-RESISTOR NETWORK 

To have exact results and to test the reliability of the 
EMT, we solve numerically a system of resistors on a fi- 
nite square lattice. Each link — (i + l,j) along x 
or — + 1) along y of the lattice is character- 
ized by a resistance which can vanish at a local critical 
superconducting temperature T c . A given fixed voltage 
is applied at two opposite edges of the square N x N 
cluster, while the two other edges have open boundary 
conditions. Given a random distribution of resistivities 
and/or T c 's on the links, the system adjusts the voltage 
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at each node and the currents along each link by im- 
plementing Ohm's law on the links and Kirchoff's law for 
current conservation on the nodes. These laws provide a 
set of 3N 2 —2N linear equations to be solved to determine 
the iV 2 — 2N voltages of the nodes (the voltages on two 
sides of the cluster are fixed) and the currents of the 2N 2 
links (the inside N 2 — 2N links plus the 2N incoming and 
outgoing links). Of course the system could be reduced 
to a smaller number of equations by simply inserting by 
hand the expressions of the voltages or of the currents 
obtained by Ohm's or Kirchoff's laws. However, in order 
to get a more transparent mapping of currents and po- 
tentials on each link or node, we choose this straightfor- 
ward representation. We proceed as follows. We extract 
from a given distribution (which can be spatially uncor- 
related or correlated) a critical temperature on each link. 
Given T c on a link, the normal-state resistivity is deter- 
mined. It can be taken as independent from T c , or it 
can be larger for less metallic regions (i.e., those having 
a smaller T c ). The latter case is aimed to account for 
the effect of microscopic impurities on the critical tem- 
perature, and it will be discussed in Sec. IV-B within 
the context of the Finkelstein theory^ The resistors on 
the links can even display a temperature dependence to 
describe, e.g., the contribution of superconducting fluc- 
tuations to the reduction of R upon approaching the local 
T c . This case will be discussed in Sec. IV-A within the 
context of both Gaussian (Aslamazov-Larkin) or vortical 
(Berezinski-Kosterlitz-Thouless) superconducting fluctu- 
ations. In any case, at a given temperature a given 
random realization of T c 's gives rise to a realization of 
normal-state resistivities. Then the set of linear Ohm's 
and Kirchoff's equations is numerically solved for a Nx N 
clusters with N ranging from 50 to 200. The total cur- 
rent / flowing at one edge of the cluster is evaluated by 
summing the currents of the N horizontal links (if the 
overall constant voltage difference V is applied to the 
two equipotential sides at x = 1 and x — N). The ra- 
tio V/I determines the total resistance of the cluster at 
any given temperature for the specific random realization 
of T c 's and of the related local resistances in each bond 
R{T C , T). This exact numerical solution is then compared 
with the solution of the EMT. 



IV. SPATIALLY UNCORRELATED 
DISTRIBUTIONS OF CRITICAL 
TEMPERATURES 

We start considering the case of spatially uncorre- 
cted distributions of critical temperatures. As we men- 
tioned above, we will consider a bimodal distribution as 
a possible realization of smoothly vanishing resistivity 
curves. The T c distribution is then formed by two gaus- 
sian peaks: the first one with weight w\ is centered at 
a positive T c = Tj representing metallic regions with 
finite superconducting temperatures, while the second 
Gaussian with a weight W2 represents regions with non- 





AT/T C 


1 


0.18 


0.75 


0.28 


0.5 


0.5 


0.4 


0.71 



TABLE II. Transition width for the bimodal distribution, ac- 
cording to the results shown in Fig. [2] 



superconducting character and it is located at negative 
T c 's (the precise location of this second Gaussian is im- 
material as long as its tail is negligible on the positive T c 
side). Typically we chose T\ = 1 and o\ s» 0.14 — 0.16. 
The relative weights W\ and W2 = 1 — W\ of the two 
components of the bimodal distribution tunes the ratio 
between regions which can and cannot become supercon- 
ducting. Of course the case of a simple Gaussian distri- 
bution of superconducting temperatures is recovered by 
choosing w\ = 1. Fig. [2] reports the case where all the re- 
gions (both superconducting and non-superconducting) 
have the same normal-state resistivity po = 1. As soon as 
T decreases below the local T c , the local resistor becomes 
superconducting and the local resistance is switched off. 
Starting from the purely Gaussian case w\ = 1 (black 
solid lines) we progressively reduce the weight of the 
metallic-superconducting regions to Wi = 0.75 (online 
red lines), w\ = 0.5 (online green lines), and w\ = 0.4 
(online blue lines). This has two effects on the p(T) 
curves: (i) the critical temperature decreases, and (ii) 
the superconducting transition width increases, and the 
p(T) displays a finite tail which eventually saturates to 
a finite value when w\ < 0.5. By defining the transi- 
tion width as done in Fig. [T] one obtains the values re- 
ported in Table II. When w\ = 1 the transition width can 
be estimated from Eq. ([3]), that gives the slope at the 
transition p' em (T c ) = 2p Q W(T c ) = y/2/npo/ai, so that 
AT C /T C — y/n/2(ai/Ti). Since we are using a distribu- 
tion having ax/Tx = 0.14 we obtain AT C /T C « 0.18, as 
reported in Table II. However, when wi decreases AT/T C 
increases, and simultaneously p(T) displays a positive up- 
ward curvature after the regime of linear decrease. As far 
as the comparison between the EMT (the dashed lines) 
and the numerical exact results is concerned, as it is ap- 
parent in Fig. [2 they coincide at any value of Wi, except 
in a very small critical region around the transition. This 
is quite remarkable since it is commonly believed that 
EMT works well only for dilute systems or for binary 
alloys of chemical species with similar resistivity^ Here 
the situation involves, instead, a binary alloy with finite 
(= po) and vanishing resistances at any relative concen- 
tration. We next explore the performance of EMT relax- 
ing the assumption of constant normal-state po for two 
cases of physical interest. In the first case a specific tem- 
perature dependence of the local resistances arises from 
paraconductive superconducting fluctuations inside each 
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FIG. 2. (Color online) Resistance curves for bimodal gaussian 
distribution of T c 's. The solid lines are the exact results, the 
dashed ones report the EMT calculation. The upper scale is 
the deviation of T from the average T c in units of <n. The 
various colors refer to different weights of the two peaks of 
the bimodal distribution: black w\ = 1; red wi = 0.75; green 
wi = 0.5; blue Wi = 0.4. The upper peak of the bimodal 
distribution has an average T c = 1 and variance a — 0.14. 



grain. In the second case quenched disorder induces a re- 
lation between the normal-state resistivities and the local 
critical temperature. 



A. Resistance with paraconductivity fluctuations 

1. Gaussian superconducting fluctuations 

In the case of gaussian superconducting fluctuations 
the well-known theory by Aslamazov and Larking in- 
troduces a power-law correction to the normal-state con- 
ductivity arising from Cooper-pair fluctuations. This in- 
crease of the conductivity diverges at T c with a pref- 
actor which is universal in two dimensions even in the 
presence of disorder. Although other corrections to the 
conductivity become influential in the case of s-wave su- 
perconducting fluctuations, namely the Maky-Thompson 
contribution and the density of state corrections^ in 
two dimensions the overall importance of the paracon- 
ductive corrections is well estimated by the universal 
Aslamazov-Larkin coefficient. This gives rather precise 
indications on the strength of the paraconductivity cor- 
rections and allows to avoid unrealistic assumptions on 
the role of paraconductive fluctuations in shaping the re- 
sistivity curves. Fig. [3] reports two cases in which the 
normal-state resistivity is lowered by paraconductivity 
effects according to the Aslamazov-Larkin theory of su- 
perconducting fluctuations. In this case 



when the local T c is lower than T, while it is zero as soon 
as T becomes lower than T c . Here po(T) is the ratio of the 
resistivity to its normal-state value (so that po (T) = 1 at 
high T), and within the ordinary AL theory the param- 
eter a is given in terms of the normal-state sheet resis- 
tance by a = Rn/Rc- Having in mind the experiments 
in the superconducting interfaces, where normal-state re- 
sistivity ranges bewteen 0.5 and 2 kf2^r— we use as an 
indicative value a = 0.05. In Fig. [3ja) the critical tem- 
peratures have a Gaussian distribution, while Fig. (3][b) 
reports the case of a bimodal distribution with weight 
W\ = 0.5 of the gaussian peak with finite average T c . 
Again the EMT (dashed lines) reproduces well the ex- 
act numerical data (symbols) in both cases. We notice 
in passing that, although the paraconductivity correction 
diverges in each grain (i.e., in each resistor approaching 
its T c ), the effect on the whole system is vanishingly small 
around the global T c . This can be understood thinking 
that around the global T c only a minor fraction of the re- 
sistors is still normal: for the whole system it makes little 
difference whether the very last few resistors of the per- 
colating cluster have their resistance constant or lowered 
by AL fluctuations. In general we notice that the para- 
conductivity corrections do lower the slope of the resistiv- 
ity curves thereby increasing the width of the global su- 
perconducting transition (compare squares and circles). 
However, it is also clear that assuming a realistic para- 
conductivity strength one cannot account for the experi- 
mentally observed transition widths without considering 
also a substantial effect of disorder (see diamonds). Re- 
sistivity slopes as small as the observed ones (see Table 
I) can only occur if the width is ruled by disorder (i.e., 
by the width of the T c distribution), while AL paracon- 
ductivity can at most provide a moderate increase of the 
width. 



2. Vortical superconducting fluctuations 

Owing to the two-dimensional character of the super- 
conducting films, we also consider the possibility that 
paraconductivity corrections might arise from vortical 
fluctuations^ While the functional form of these cor- 
rections in terms of the coherence length £ is the same 
as for the Gaussian fluctuations, i.e., 5a oc £ 2 , the 
prefactor is not universal and the temperature depen- 
dence reflects the exponential behavior of the coher- 
ence length in the Berezinski-Kosterlitz-Thouless (BKT) 
transition ; 23 ! 24 i.e., 
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where t = (T — Tbkt) /Tbkt is the distance from the 
BKT transition temperature Tbkt, A is a constant of 
order one, and the parameter b is controlled by the dis- 
tance bewteen the Tbkt and the mean-field temperature 
T®, above which Gaussian fluctuations are restored. As it 






FIG. 3. (Color online) Resistance curves in the presence of 
AL paraconductivity fluctuations and/or disorder. Symbols 
represent the exact results, the dashed lines report the EMT 
calculation, (a) Gaussian distribution of disorder. (Squares) 
Gaussian distribution of T c 's with a = 0.16 and T c = 1 and 
no AL fluctuations (a — 0); (diamonds) resistance in the pres- 
ence of AL fluctuations (a = 0.05), but no disorder [a = 0); 
(circles) resistance in the presence of both AL fluctuations 
(a — 0.05) and disorder (a — 0.16 and T c = 1). (b) Same as 
in (a), but with a bimodal gaussian distribution of T c 's with 
Wi = 0.5. 

has been discussed in detail in Ref. [24| b = 2ay/t^, where 
t c = (T° — Tbkt) /Tbkt and a is a measure of the en- 
ergy of the vortex core, expressed in units of the value it 
assumes within the standard JTY-model approach to the 
BKT transition. Having in mind these definitions, we 
first explored the effect of local vortical paraconductiv- 
ity fluctuations on the global superconducting transition 
using t c = 0.1, a = 0.2 (i.e., b = 0.126), and A = 2, 
which correspond to realistic parameter values, accord- 
ing to the estimates for thin disordered superconducting 
films in Ref. [25|. The results are shown in Fig. |4] where, 
similarly to what found above for the Gaussian paracon- 
ductivity effects, we find that the most realistic estimate 
of the parameters produces only rather small corrections 
to the widths of the transition. It is quite apparent from 
Figs. [3] and E] that, for a realistic choice of parameters, 
the effects of AL fluctuations and vortical fluctuations 
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FIG. 4. (Color online) Resistance curves in the presence of 
vortical paraconductivity fluctuations and/or disorder. The 
symbols represent the exact results, the dashed lines report 
the EMT calculation, (a) Gaussian distribution of disorder. 
(Squares) Gaussian distribution of T c 's with a — 0.16 and 
no vortical fluctuations (A = 0); (diamonds) resistance in the 
presence of vortical fluctuations [t c — 0.1, a = 0.2 (6 = 0.126), 
and A = 2], but no disorder (a = 0); (circles) resistance in 
the presence both of vortical fluctuations [t c = 0.1, a = 0.2 
(b = 0.126), and A = 2] and of disorder (a = 0.16). (b) Same 
as in (a), but with a bimodal gaussian distribution of T c 's 
with iui = 0.5. 



are quite similar in the absence of disorder as well as in 
the presence of both a gaussian or a bimodal distribution 
of T c 's. 

Since in the literature the correct estimate of the pa- 
rameters within the Halperin-Nelson formula (|6]) is of- 
ten debatedj2i we decided to explore further the limit 
of the BKT theory in reproducing correctly experimen- 
tal data in superconducting interfaces. In particular, we 
considered also a set of parameters able to yield a larger 
transition width, regardless of their microscopic deter- 
mination. In Fig. [5] we show the results for two set of 
Halperin-Nelson parameters and the related fit to data 
from an unbiased sample in Ref. 0. In panel (a) we re- 
port the results for a set (HN1) with t c = 0.1, a = 0.4 
(b = 0.25), and A = 2, in the presence of a gaussian 
distribution of T c 's (wi — 1) with o\ — 0.035 K and 
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FIG. 5. (Color online) (squares) Resistivity data of a 
LaA103/SrTi03 sample at zero bias from Ref. 0. (a) (Black) 
Solid line, fit in the presence of both gaussian disorder (wi — 
1, c i = 0.035 K, and T c = 0.195 K) and vortical fluctuations 
(HN1) with t c = 0.1], a = 0.4 (6 = 0.25), and A = 2. The 
(blue online) dashed line is for a similar parameter set (HNf), 
with t c = 0.1, a = 0.4 (b = 0.25), and A = L8, and a bimodal 
distribution Wi = 0.5, cti = 0.025 K, and T c = 0.23 K. (b) 
Squares, same as in (a); (black) solid line, fit in the presence 
of vortical fluctuations only (HN2) with t c = 0.7, a — 0.8 
{b = 1.34), A = 0.71, and T c = 0.157 K. 



T c = 0.195 K (black solid line) and for a similar set 
(HNT) with t c = 0.1, a = 0.4 (b = 0.25), and A = 1.8, in 
the presence of a bimodal distribution of T c 's (u>i = 0.5) 
with cri = 0.025 K and T c = 0.23 K (blue dashed line). 
In panel (b), instead, we attempt to fit the same data 
with vortical fluctuations only (i.e., a = 0). Of course we 
need a quite more substantial amount of vortical fluctua- 
tions as given by a different (rather unrealistic^-.) param- 
eter set (HN2), with t c = 0.7, a = 0.8 (b = 1.34), and 
A = 0.71. It is apparent that vortical fluctuations alone 
can account for broad transitions [and also produce the 
tail characteristic of the LaA10 3 /SrTi0 3 (LAO/STO) or 
LaTi0 3 /SrTi0 3 (LTO/STO) interfaces]. However, the 
parameters needed to produce the fit (which, by the way 
is quite poor in reproducing the high-temperature down- 
ward curvature) can hardly be justified on microscopic 
grounds. On the other hand, a moderate amount of dis- 
order (notice that T c « 0.2 K is in physical units here 
and therefore ct\ 0.03 K is comparable with the values 
<ti/T c 0.15 of the figures, having normalized T c = 1) 
allows quite reasonable fits using a far more standard 
set of Halperin-Nelson parameters [see Fig. [Sl[a)]. We 
find also in this case that the overall width of the transi- 
tion can be captured by the gaussian distribution of T c 's 
(with a moderate contribution of superconducting fluc- 
tuations), but to reproduce the low-temperature tail one 
needs a bimodal distribution. 



B. Resistance with quenched disorder 
(Finkelstein's theory) 



Another physically interesting case occurs when 
quenched impurities affect the critical temperature of the 
grains. In Ref. [9j a connection was established between 
disorder and the critical temperature in two-dimensional 
samples by deriving a relation between the resistance and 
the critical temperature in the presence of both disorder 
and Coulomb repulsion. Here, assuming that each grain 
is large enough to make the theory of Ref. applica- 
ble, we extract randomly (from a Gaussian distribution) 
the local critical temperatures of the resistors and we 
consequently assign the local resistivity. In practice we 
numerically invert the relation 



Tc 
TP 



7 - t/A + ^Jt/2 



7-£/4- v/t/2 



(7) 



where, following the notation of Ref. |9j, T® is the criti- 
cal temperature of the clean system, t = Rn/R' c , 7 = 
1/Iti(T®t), t is the elastic scattering time, and the re- 
sistance R' c = 2n 2 h/e 2 = 81.1 Kf2. For the sake of 
concreteness, in our calculation we use 7 ~ —0.11, ob- 
tained choosing T c ° = 230 mK as the superconducting 
temperature for the clean system and r = 5 x 10~ 7 mK -1 
for the value of the elastic scattering time. These val- 
ues have been chosen to reproduce the correct orders of 
magnitude of experiments in LTO/STO samples^ The 
resulting T c vs. The Rn curve is reported in the inset 
of Fig. [6] The results reported in Fig. [6] again display a 
quite good agreement between the exact numerical cal- 
culations and the EMT calculations (the barely visible 
dashed lines). We notice here that the progressive in- 
crease of the resistivity upon considering more disordered 
systems with lower T c 's is associated with an increasing 
slope of the resistivity curves near T c . By rescaling all 
the curves to let them assume the same high-temperature 
value pa(T 3> T c ) = 1, the various curves are found to 
acquire quite similar slopes. This intriguing feature of 
our data is present in the case of LTO/STO interfaces, 
where, however, tails at T ~ T c are also found, which are 
missing in our calculations of Fig. [BJ carried out with a 
single gaussian distribution of local T c 's. 



V. SPATIALLY CORRELATED 
DISTRIBUTIONS OF CRITICAL 
TEMPERATURES 



In all above cases the distribution of T c 's and the con- 
sequent distribution of the resistances were spatially un- 
corrected. This generically allows for a remarkably good 
performance of the EMT. Here, instead, we challenge 
EMT in the case of various space correlations. 
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FIG. 6. (Color online) Resistance for systems with gaussian 
distribution of local T c 's centered at different average T c 's 
and related local resistances from Eq. (0 with R c = 81.1 Kf2, 
T c ° = 230 mK, and r = 5 x 1CT 7 mK -1 (the inset reports 
the T c vs. p curve corresponding to Eq. [7J. The variance 
of the gaussian distribution is always a = 30 mK. Symbols 
refer to exact calculations while dashed lines report the EMT 
curves. The gaussian distributions are centered at T c = mK 
(blue circles and line), T c = 70 mK (green squares and line), 
T c — 150 mK (orange diamonds and line), and T c = 200 mK, 
(red triangles and line). 



A. Short-range averaged critical temperatures 

In order to investigate the effects of short-range space 
correlations, we implement a short-range averaging pro- 
cedure of the critical temperatures. As a first step, we 
extract the local T c 's of the bonds on our cluster from a 
given gaussian distribution. Then, the value of T c on each 
bond is averaged with the values on the six nearest neigh- 
bor bonds. Of course, this yields gaussian distributed 
T c 's with a variance reduced by a factor yfl and, at the 
same time, creates space correlations on distances of the 
order of two lattice spacings. This short-range correla- 
tion can be extended by iterating the averaging proce- 
dure: the once- averaged T c 's can be averaged again with 
the six nearest neighbors, leading to space correlations 
up to four lattice spacings. The results we display below 
are always obtained with a two-step averaging protocol. 
Calculations with one- or three-step averaging procedures 
yield similar results. Fig. [7J reports with the solid black 
dots a calculation obtained starting from a Gaussian dis- 
tribution with variance a = 0.33, centered at T c = 1. 
The corresponding result for a EMT calculation is given 
by the black dashed line. Since EMT completely ignores 
the space correlations, the results are the same as those 
obtained with the starting uncorrelated set of random 
T c 's, but for a trivial rescaling of the variance. Here, we 
numerically find that the variance entering the EMT ex- 
pression [cf. Eq. (@|] with w 2 = is tj\ — 0.09. Clearly 
the averaging procedure spatially correlates the high-T c 



FIG. 7. (Color online) Resistance of a short-range correlated 
random resistor network with T c — 1 and a starting a = 0.33. 
The solids dots are the exact results, while the dashed line is 
obtained from numerical evaluation of the EMT, Eq. [1] with 
a = 1. The red dot-dashed line is instead obtained from Eq. 
H with wi = 1, W2 = 0, a = 1.6, Ti = 0.997, and ai = 0.09. 



regions, which form short-range clusters. The inset in 
Fig. [7] displays a 100 x 100 cluster in which the black re- 
gions correspond to regions where the local T c is larger 
than the average value T c = 1. This lowers the exact 
resistance curve with respect to the EMT one, since the 
superconducting cluster is formed by objects with an ef- 
fectively higher connectivity. Of course, this is true only 
far from percolation, because on a large scale it is instead 
immaterial whether the percolating objects are the sin- 
gle bonds or these small short-range clusters. Moreover, 
it is clear that the averaging procedure does not break 
the symmetry between high-T c and low-T c bonds (even 
if the variance is made smaller). Therefore, the perco- 
lation threshold stays the same and it is reached in two 
dimension when half of the bonds are superconducting. 
As a consequence, the curve of the exact calculation fol- 
lows at high temperature the behavior of the EMT with a 
higher connectivity [this is represented by the (red online) 
dot-dashed line obtained from Eq. (Q}, with an effective 
connectivity numerically estimated to be a ~ 1.6]. On 
the other hand, by lowering the temperature, the resis- 
tance is dominated by percolation effects, and a change 
of curvature occurs in p(T), which vanishes at T = 1, 
corresponding to a concentration 1/2 of superconducting 
bonds. This effect clearly introduces a "tailish" charac- 
ter in the low-T part of the resistivity, which, however, 
is not large enough to reproduce quantitatively the data 
in LAO/STO and LTO/STO interfaces. 



B. Patches 

Besides averaging the T c 's over the first neighbors, we 
also introduce space correlations by simulating the occur- 
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FIG. 8. Resistance curves for a 100 x 100 system with patchy 
structure (see lower inset). Each of the 80 patches (of radius 
r — 5 lattice units) has a uniform local T c extracted from a 
Gaussian distribution with T c = 1 and a = 0.16. The un- 
derlying matrix has local T c 's extracted from a similar Gaus- 
sian centered at T c = —100. (negative T c 's are cumulated in 
T c = 0). The solid line is the exact result, while the dashed 
line is the EMT calculation. Upper inset: Total distribution 
of T c 's. 



rence of patches, where the T c 's are similar in a region 
of given size. We first select a set of sites randomly dis- 
tributed in our N x N cluster. Around each of these 
"seed" sites we define a patch of a given radius r. Then, 
we extract the critical temperatures inside the patches 
from a given gaussian distribution centered at a higher 
value of T c , while the bonds outside the patches have 
vanishing T c . In this way, we aim to simulate a bimodal 
distribution with spatial correlation over a range r typical 
of a sample where different extended regions have a more 
or less marked superconducting character. The lower in- 
set of Fig. [8] represents the case in which 80 patches with 
r = 5 are introduced in a 100 x 100 cluster. The main 
panel of Fig. |8] reports the related resistances both from 
exact (solid line) and EMT (dashed line) calculations. In 
both cases the system displays an initial rapid decrease of 
p, which then saturates at finite values upon decreasing 
T. While the EMT and exact data are quite similar at 
high temperatures, the low-T saturation values are dif- 
ferent. In particular the prediction of EMT is smaller 
than the exact result. Indeed, the exact solution is more 
resistive because the superconducting bonds are grouped 
into patches thereby leaving more extended surrounding 
regions, which are resistive. As a limiting case, one could 
imagine to form a single patch with weight 1/2 located 
inside the square cluster. Since one would have 1 /2 of the 
bonds superconducting, EMT would predict the overall 
system to become superconducting. On the other hand, 
this single big patch being completely surrounded by a 
resistive region, the exact calculation would give a finite 
resistance for the system. This indicates that this real- 
ization of patches naturally results in higher resistances 
in the exact calculation, where the space correlation ef- 
fects play a role, with respect to the EMT, where only the 



FIG. 9. Same as in Fig. [8] but for a system with 180 patches. 



overall weight of superconducting bonds matters. In the 
present model, the size of superconducting regions can 
be increased either by increasing the number of patch 
centers (the "seeds") and/or by increasing the radius of 
the patches r. Fig. |9] reproduces a case with 180 patches 
with r — 5 in a 100 x 100 cluster. Clearly (see lower in- 
set), there is a large majority of highly superconducting 
sites, and the EMT predicts that the whole system be- 
comes superconducting at T c = 0.91 (black dashed line 
in the main panel). Remarkably, this prediction repro- 
duces well the exact solution given by the black solid 
line. However, close inspection of the lower inset shows 
that the superconducting character of the exact solution 
is only due to a small overlap of the patches in the cen- 
tral region of the cluster. If this small superconducting 
region were absent, the whole assembly of patches would 
not percolate and the system would display the finite re- 
sistance due to this small resistive region. The above 
patchy structures have a random character, which makes 
it difficult to unambiguously establish the effective con- 
nectivity of the correlated regions. Therefore, we also in- 
vestigated a toy model where the patches no longer form 
around random centers, but they form in a rather regular 
way. Specifically, we simply subdivided our N x N clus- 
ter into smaller (N/d) x (N/d) squares. This gives rise to 
a checkerboard support with d x d squares each having 
a given T c randomly extracted from a gaussian distri- 
bution. As a consequence, the critical temperatures of 
this coarse-grained system are still gaussian distributed 
but also display a strong space correlation, since they 
all coincide within each of the d x d square subclusters. 
We also notice in passing that to obtain a reasonably 
good Gaussian distribution quite large values of N are 
required, because the statistical distribution of T c 's is 
substantially degraded by extracting just one random T c 
inside each subcluster. We notice, however, that in this 
toy model each square subcluster is only connected to 
four other subclusters (but for a small effect in the cor- 
ners). The exact calculation (not shown here) in this 
case agrees remarkably well with the prediction of the 
EMT. This demonstrates that connectivity plays a much 
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more important role than the simple space correlation. 
Of course, whenever this latter influences the connectiv- 
ity (like in the models of Sees. V A and V B), space 
correlation again becomes a main actor of the game. 



C. Long-range spatially correlated distributions 

To complete the analysis on the reliability of EMT 
in the presence of spatial correlations, we also investi- 
gate the behavior of the exact resistance and EMT when 
long range correlations are present. We empirically sim- 
ulate this situation by adjusting the seeds of the random- 
number generator producing the distribution of random 
T c 's, while visiting the NxN lattice sites, so that the T c 's 
are highly correlated along lines. In particular, one can 
easily introduce long-range correlated patterns like the 
one reported in Fig. [TUtah where high values of T c along 
diagonals alternate with low values forming a stripe-like 
texture. In this case, the exact resistance displays rapid 
drops whenever the temperature reaches a value at which 
several neighboring bonds become simultaneously super- 
conducting, substantially extending the size of the previ- 
ously formed superconducting clusters [see the thin solid 
lines in Fig. [TUTb)]. On the other hand, it then becomes 
difficult to build a fully percolating cluster and more or 
less sizable tails are generated. Lacking any information 
about the space structure of the random system, EMT 
fails in reproducing this situation. This is clearly visi- 
ble comparing the two blue lines in Fig. llOf b). reporting 
the exact (solid) and EMT (dashed) results for the same 
specific random realization of Fig. HUT a). 

The sharp drops occurring in the individual disordered 
realizations can however be smoothed by averaging over 
various realizations of the long-range correlated patterns 
(simulating, for instance, a system with differently ori- 
ented textures). The thick solid (red online) line in Fig. 
UOf b) is indeed obtained by averaging over the thin solid 
resistance curves obtained from seven different realiza- 
tions with differently oriented and correlated diagonal 
stripes like the one reported in Fig. HUT a) . As mentioned 
above, each of these realizations has a resistance quite 
different from that obtained with EMT. The (blue on- 
line) dashed line is the the EMT result obtained with 
the same gaussian distribution of local T c 's used in Fig. 
EIKwi = 1, w 2 = 0, Ti = 1, and a x = 0.1). We also fitted 
the average curve [see inset in Fig.fTUTb)] by the analytic 
form in Eq. ((4]). This time the fit indicates an effective 
reduction of the connectivity, with a — (i.e., D = 1), 
T c = Ti = 1.15, cri = 0.075. While it is rather natural 
that some generic one-dimensional character arises from 
spatial correlations like in Fig. [TUT a). the fact that the 
averaging procedure transforms the various quite differ- 
ent resistances of the individual realizations into a single 
nearly one-dimensional resistance curve is interesting. 



(a) + 
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FIG. 10. (Color online) (a) Example of a highly spatially 
correlated realization of random Gaussian distribution of T c 's 
with T c = 1 and a = 0.1. (b) Resistance curves for various 
realization of similar random spatially correlated T c 's (thin 
solid black lines) and their average (thick solid red line). One 
of these specific curves has been drawn with a blue dot-dashed 
line together with its numerically calculated EMT (dashed 
line). Inset: The average resistance and a fit from Eq. Q 
with wi = 1, W2 = 0, a = (corresponding to an effective 
dimensionality D = 1), Ti = 1.155, and <ti = 0.075. 



VI. CONCLUSIONS 

In this paper we aimed to understand whether meso- 
scopic inhomogeneities could provide the physical mecha- 
nisms producing unusually broad superconducting tran- 
sitions in various nearly two-dimensional systems, like 
the quasi-two-dimensional electron gases formed at the 
interfaces between SrTiOs and LaA103 or LsTiOa. Fur- 
thermore we wanted to get insights on the "tails" appear- 
ing in the resistance curves of LTO/STO and LAO/STO 
interfaces. By modeling the inhomogeneous system by a 
random resistor network, we systematically investigated 
various realizations of random resistors, which we solved 
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both by an exact numerical approach and by a standard 
EMT approach. We started with the simple assump- 
tion of a random distribution of local T c 's with global su- 
perconductivity occurring when superconducting regions 
(links) will eventually percolate by lowering the temper- 
ature. We then discussed various physical ingredients to 
account for the main experimental features of the super- 
conducting transition at the interfaces. First of all, the 
transitions in these systems are anomalously broad and 
we clarified that this can hardly be attributed to stan- 
dard gaussian or vortical (i.e., BKT-like) paraconduc- 
tivity fluctuations. This is particularly evident for the 
AL paraconductivity, but it also holds for the Halperin- 
Nelson formula. Indeed, the set of parameters to be used 
to fit the experiments without disorder comes out to be 
rather unphysical. Therefore, superconducting fluctua- 
tion effects can contribute to, but cannot completely ac- 
count for, the substantial width of the transition and 
for the small slopes of R(T) around T c . On the con- 
trary, within our model the width of the transition mainly 
stems from the width of the distribution of T c 's among 
the various superconducting mesoscopic regions. It turns 
out that a typical gaussian distribution of local T c 's has 
the right downturn to mimic the experimental resistance 
curves, which display a moderate rounding above T c fol- 
lowed by a rather broad, nearly linear decrease when T 
is lowered. 

We also found in passing that the degree of reliability 
of EMT in solving different disorder realizations is high 
and standard EMT generically reproduces quite well the 
exact resistance curves, whenever the space correlations 
between the mesoscopic regions are negligible. This re- 
sult is general, and is not limited to dilute systems or 
to mixtures of similar resistances. We rather find that 
EMT works well for mixtures of metallic (finite R) and 
superconducting (i.e., R — 0) resistances, for any filling 
ratio, possibly including the effects of superconducting 
fluctuations. On the contrary, we find that increasing 
the space correlations leads to more substantial failures 
of the EMT. 

An additional interesting and specific feature of the 
two-dimensional electron gas formed at the interface of 
the SrTi03 substrate is a more or less long tail in the 
low-T part of the R(T) resistance curves. These tails 
can even flatten to form a finite- R plateau in cases where 
the system stays non superconducting despite a substan- 
tial decrease of R(T) indicating the presence of a siz- 
able superconducting fraction of the two-dimensional gas. 
We found that this feature is not easily reproduced and 
we only found few specific cases, where it could be re- 
produced within our model. In the case of non spa- 
tially correlated local T c 's we found that a tail is only 
present when a bimodal distribution of T c 's is assumed, 
with (or very near to) the specific 1/2 — 1/2 distribu- 



tion of the relative weight of the two distribution peaks. 
Since in the real systems the tails are observed over a 
broad range of fillings, magnetic fields, bias, it is hard to 
understand how these different conditions can all corre- 
spond to such specific 1/2 — 1/2 constraint on the dis- 
tribution of the mesoscopic local T c 's. The 1/2 — 1/2 
condition can be rephrased by saying that at perco- 
lation, the whole set of superconducting domains ran- 
domly distributed in the two-dimensional system has to 
be exhausted^ One possibility which deserves to be ex- 
plored is that some microscopic mechanism acts to create 
some low-dimensional skeleton on which superconductiv- 
ity takes place. Localization effects have already been 
proposed as a possible mechanism to generate a frac- 
tal substrate for superconductivity^ and a glassy super- 
conducting phasei 15 ' 16 If superconductivity occurred on 
this low-dimensional substrate, most if not all the bonds 
of the substrate have to be superconducting before the 
whole "skeleton" acts as a whole superconducting path. 

A further mechanism to explore in order to reproduce 
tails in the resistance of the quasi-two-dimensional elec- 
tron gas is provided by space correlations (indeed, the 
above scenario of fractal or glassy superconducting phase 
is a form of space correlation) . When space correlations 
are short-ranged (see Sec. V A), we detect the presence 
of a short tail (see Fig. [7]). One can expect that extend- 
ing the strength and the range of these correlation could 
result into a strengthening of the tail feature. This ex- 
pectation is confirmed by investigating space structures 
where more metallic domains cluster forming patchy tex- 
tures embedded into less metallic matrices with very low 
or zero superconducting temperatures (see Figs.[8]and[9]). 
In these cases again bimodal T c distributions are formed 
and again "tailish" resistances are obtained whenever the 
weight of the high-T c distribution peak approaches one 
half (or less when a plateau is obtained). 

Finally, to reproduce tailish resistance curves, we ex- 
plored the case of textured domains (forming, e.g., stripe- 
like paths) with long-range space correlations. If the real 
system is formed by a patchwork of randomly oriented 
stripe domains, one would obtain resistance curves with 
substantial tails (see Fig. [TU)) . Indications of a spatially 
structured system hosting the two-dimensional electron 
gas have indeed been obtained from scanning tunneling 
microscopy experiments on LAO/STO samples^ sug- 
gesting that also this possibility is worth being further 
explored. 
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